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The brick wall model is a semi-classical approach to understand the microscopic origin of black hole 
entropy. In this approach, the black hole geometry is assumed to be a fixed classical background on 
which matter fields propagate, and the entropy of black holes supposedly arises due to the canonical 
entropy of matter fields outside the black hole event horizon, evaluated at the Hawking temperature. 
Apart from certain lower dimensional cases, the density of states of the matter fields around black 
holes cannot be evaluated exactly. As a result, often, in the brick wall model, the density of states 
and the resulting canonical entropy of the matter fields are evaluated at the leading order (in terms 
of h) in the WKB approximation. The success of the approach is refiected by the fact that the 
Bekenstein-Hawking area law — viz. that the entropy of black holes is equal to one-quarter the 
area of their event horizon, say, Ayi — has been recovered using this model in a variety of black 
hole spacetimes. In this work, we compute the canonical entropy of a quantum scalar field around 
static and spherically symmetric black holes through the brick wall approach at the higher orders 
(in fact, up to the sixth order in K) in the WKB approximation. We explicitly show that the 
brick wall model generally predicts corrections to the Bekenstein-Hawking entropy in all spacetime 
dimensions. In four dimensions, we find that the corrections to the Bekenstein-Hawking entropy 
are of the form [^JJ log^j^] , while, in six dimensions, the corrections behave as 4- A^^ logVljj] , 
where {m,n) < 1. We compare our results with the corrections to the Bekenstein-Hawking entropy 
that have been obtained through the other approaches in the literature, and discuss the implications. 
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I. INTRODUCTION AND MOTIVATION 

The concept of black hole entropy was originally intro- 
duced by Bekenstein [1~4] to resolve certain thcrmody- 
namical paradoxes that arise in the presence of the black 
holes and, in particular, to preserve the universal applica- 
bility of the second law of thermodynamics. Soon after 
Bekenstein's proposal — based on their classical, macro- 
scopic behavior — the thermodynamic properties of black 
holes were formalized as the four laws of black hole me- 
chanics [5]. Specifically, it was argued that, as the area 
theorem of classical general relativity closely resembles 
the statement of the second law of thermodynamics, the 
area of the black hole event horizon {A-^) can be inter- 
preted as the physical entropy associated with the black 
hole. This association, in turn, led to the identification 
of the surface gravity (k) of the black hole (which, for a 
stationary black hole, is a constant all over the horizon) 
as the temperature of the hole. 

The laws of black hole mechanics were placed on a firm 
footing when, a year or two later. Hawking [6, 7] showed 
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that, in the presence of quantum matter fields, a body 
that collapses into a black hole emits thermal radiation 
at the temperature 




where h, c and denote the Planck constant, the speed 
of light and the Boltzmann constant, respectively. The 
above Hawking temperature fixes the constant of pro- 
portionality between the temperature of the black hole 
and its surface gravity and, therefore, between the en- 
tropy and the area of the hole. One finds that the en- 
tropy of black holes are given by the following Bekenstein- 
Hawking area law: 




r /2 

where £-p^ = {Gh/c'^ denotes the Planck length with 
G being the Newton's constant. 

Black hole entropy assumes considerable importance 
due to the fact that it may provide us with an in- 
sight to the microscopic structure of the gravitational 
theory through the microcanonical, Boltzmann relation 
S — (/Cg lnf2), where is the total number of quantum 
states that are accessible to a black hole that is described 
by a small set of classical parameters. The different 
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approaches that have been adopted in the literature to 
understand the microscopic origin of black hole entropy 
can be broadly classified into two categories, (i) Count 
the "microstates" by assuming a fundamental structure 
like D-branes, spin networks or conformal symmetry [8- 
11]. (ii) Associate the black hole entropy to the quan- 
tum fields propagating in the fixed black hole spacetime. 
Count the microstates of these quantum fields [12-17]. 

Although none of the above approaches can be consid- 
ered to be complete; all of them — within their domains 
of applicability — by counting certain microscopic states 
yield the semiclassical result (2) in all spacetime dimen- 
sions d > 3. However, all these approaches seem to lead 
to different sub-leading corrections to the Bekenstein- 
Hawking entropy. For instance, (i) the prefactor to the 
logarithmic corrections obtained using the spin-networks 
and conformal symmetry [18-22] are different from the 
one obtained using the statistical fluctuations around 
thermal equilibrium [23]. (ii) the power-law corrections 
obtained using the Noether charge approach [15] are dif- 
ferent from those via entanglement of the modes between 
inside and outside the horizon [24]. In other words, even 
though different degrees of freedom lead to the universal 
Bekenstein-Hawking entropy — quite naturally — they 
lead to different sub- leading terms. This indicates that 
the key to the understanding of the statistical mechanical 
interpretation of Bekenstein-Hawking entropy may lie in 
the origin of the sub- leading contributions. 

Physically, it is natural to expect corrections to (2): 
Bekenstein-Hawking entropy is a semi-classical result and 
there are strong indications that this is valid for large 
black holes [i.e. when horizon radius is much larger than 
Planck length]. However, it is not clear, whether this 
relation will continue to hold for the Planck size black 
holes. Besides, there is no reason to expect that the 
Bekenstein-Hawking entropy to be the whole answer for 
a correct theory of quantum gravity. 

In this work, we calculate the higher order WKB con- 
tributions to the Bekenstein-Hawking entropy from the 
brick- wall model [12, 25, 26]. We extend the zeroth-order 
(hP) WKB analysis to higher order and show that (i) The 
contribution to the entropy from the higher-order WKB 
modes is of the same order as the leading order WKB 
modes. In other words, our analysis shows that it may 
be incomplete to calculate the contribution only from the 
leading order WKB modes, (ii) The brick-wall entropy 
{S^^) leads to generic corrections to area of the form: 

5bw = -^BH + G{An) + HAh) log , (3) 

where ^(^h), and J^{Ab_) are polynomial functions of 
A^. In the case of four-dimensions, we show explicitly 
that the brick-wall entropy (upto sixth-order) has the 
form given above with Q{Aii) — 0. In the case of six- 
dimensions, G{Ah) 7^ 0. (iii) We show that, only in the 
case of Schwarzschild, J^{Ah) is a constant. 

The brick wall approach is a semi-classical approach, 
wherein the background geometry is assumed to be a 



fixed classical background in which quantum fields prop- 
agate. The entropy of the black hole is identified with 
the statistical mechanical entropy arising from a thermal 
bath of quantum fields propagating outside the horizon. 
The entropy computed in this way turns out to be pro- 
portional to the area of the horizon. This approach has 
been very popular in obtaining the leading order to the 
black hole entropy in different dimensions (for an incom- 
plete list of references, see Rcfs. [27- 42]). 

The problem with the brick- wall model [as is the case 
of any semiclassical approach] is that due to the infinite 
growth of density of states close to the horizon, one has 
to impose ultra-violet cutoff near the horizon and hence, 
the brick-wall entropy depends on the cut-off scale. [See 
section (III), for discussion on various aspects of brick- 
wall model.] Clearly, this is an undesirable feature. How- 
ever, there are several advantages of the brick- wall model 
over other approaches: (i) Unlike the Noether charge ap- 
proach [15], the brick- wall entropy depends only on the 
kinematical properties of the metric close to the horizon 
and does not depend on the dynamics. Hence, the brick- 
wall entropy (and the corrections computed in this work) 
can directly be mapped to the horizon properties. In the 
case of Noether charge approach, since, such a mapping 
is not possible, the power-law corrections does not pro- 
vide any new information about the statistical mechan- 
ical properties of black hole entropy, (ii) Unlike entan- 
glement entropy [24], the brick- wall entropy can be com- 
puted analytically for any spherically symmetric space- 
times to all orders. Also, it is not possible to compute the 
entanglement entropy for spacetime dimensions c? > 4 — 
the entropy is divergent, (iii) In the conformal field the- 
ory approach [10, 19], the black hole horizon is treated 
as boundary. However, the vector fields (which generate 
the symmetries) do not have a well-defined limit at the 
horizon [44] . If one requires that vector fields generating 
symmetries be smooth at the horizon, then the central 
charge should be zero. In other words, the conformal 
field theory analysis can only be performed close to the 
horizon like a brick-wall. 

The remainder of this paper is organized as follows. 
In the following section, we shall sketch some essential 
properties of static, spherically symmetric black holes in 
arbitrary spacetime dimensions. In Section III, we shall 
discuss the assumptions and approximations involved in 
evaluating the brick wall entropy, and outline the proce- 
dure for extending the calculation to the higher orders (in 
terms of h) in the WKB approximation. In Section IV, 
in addition to the zeroth order, we shall evaluate the con- 
tributions to the brick wall entropy of four dimensional 
black holes at the second, the fourth and the sixth orders 
(in terms of H) in the WKB approximation. In Section V, 
we extend the analysis to case of black holes in six dimen- 
sions. In Section VI, we explicitly write down the results 
for a few specific black hole solutions in four and six di- 
mensions. Finally, in Section VII, after a rapid summary 
of the results we have obtained, we shall discuss as to how 
the sub-leading contributions we have evaluated compare 
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with the results obtained from the other approaches. 

Before we proceed further, let us briefly outline the 
conventions and notations we shall adopt. Wo shall, in 
general, consider a (D+2)-dimensional^, spherically sym- 
metric, black hole spacetime. We shall work with the 
metric signature (—,+,+,•••), and use the geometric 
units wherein = c = G = 1. We shall denote the 
derivative of any function with respect to the radial co- 
ordinate r of the black hole by an ovcrprimc. The quan- 
tum field $ we shall consider will be a minimally coupled 
scalar field. 



II. KEY PROPERTIES OF STATIC, 
SPHERICALLY SYMMETRIC BLACK HOLES 

Consider the following {D + 2)-dimensional static and 
spherically symmetric line element 



g{r) 

= f{r) [-dt^ + dx^] + dnl , 



(4) 

(5) 



where /(r) and g{r) are arbitrary (but, continuous and 
differentiable) functions of the radial coordinate r, dCl"^ 
is the metric on a £)-dimensional unit sphere, and 



X 



■ 



Vfir)9ir) 



(6) 



denotes the tortoise coordinate. Throughout this work, 
we shall assume that the line-element (4) contains a sin- 
gularity (say, at r = 0) and one, non-degenerate, event 
horizon (located at, say, r = r^)^. But, we shall not as- 
sume any specific form of /(r) or g{r). In the rest of this 
section, wc shall discuss some generic properties of the 
spacetime (4) near the horizon a,t r = r^. 

In almost all approaches that evaluate the entropy 
of spherically symmetric black holes, their line-element 
close to the event horizon is approximated to be that of 
a Rindler spacetime (see, for instance, Ref. [39]). For 
the line-element (4), the Rindler behavior near the hori- 
zon can be arrived at by first carrying out the following 
transformation of the radial coordinate: 



(7) 



^ A comment on this uncommon notation may be in order. We 
do not work with two time coordinates! We find it convenient 
to assume the spherically symmetric event horizon to be a D- 
dimensional sphere. The radial and the time coordinates maJce 
the dimension of spacetime to he (D + 2). 

^ Since the event horizon is a null hypcrsurfacc, its location can be 
determined by the condition (gf" dfj,N di,N) = 0, where W is a 
scalar quantity that characterizes the liypersurface. For the line- 
clement (4), A'' proves to be a function of the radial coordinate r. 
The above equation then leads to g{r^ ) = 0, and the roots of this 
algebraic equation in turn determine the location of the horizon. 



where k is a constant that denotes the surface gravity 
of the black hole and is defined as (see, for example, 
Ref. [46]) 



Igirl (fir) 
fir) 



(8) 



In terms of the coordinate 7, the line-element (4) can be 
expressed as 



ds" = -K^ 7^ di 



dj^+r'^dnl. (9) 



Close to the horizon (i.e. near r = r^), this line-element 

reduces to 



ds' 



-K^j^dt'^ + d-f^+rl dnj 



(10) 



which describes the Rindler spacetime with a horizon 
that is located at 7 = 0. It should be stressed here that 
such a behavior is exhibited by all non-degenerate black 
hole horizons in all dimensions. 

The above derivation of the Rindler line-element near 
the horizon is essentially equivalent to expanding the 
metric components /(r) and g{r) in (4) about up to 
the linear order in the Taylor series. However, we find 
that, when evaluating the contributions to the brick wall 
entropy at the higher orders in the WKB approximation, 
we need to expand the quantities /(r) and g{r) to higher 
orders as follows: 

fir) = fir J (r - rj + (^^) (r - rj^ (Ha) 

+ (^) (r--rj^ + ..., 
g{r) = g'irj (r - r J + ( ^-^) (r - rj^ (ub) 



' 9"'irJ ' 



(r-r„)3 + ... 



As we shall see, in four dimensions, in addition to the 
surface gravity of the black hole, the corrections to the 
Bckenstein-Hawking entropy Sj^^ also depend on the sec- 
ond derivative of the metric evaluated at the horizon. 
And, in six dimensions, we find that the sub-leading con- 
tributions to /Sbjj involve the third derivative of the met- 
ric as well. 

Another quantity which we shall require in our calcu- 
lations is the proper or the coordinate invariant distance 
of the brick wall from the horizon. The proper radial dis- 
tance to the brick wall, say, he, that is located at r = /i 
is given by 



r^+h 



I 



dr 



(12) 
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On using the expansion (11) for g{r) up to the second 
order in this integral, we obtain the fohowing relation 
between h and he- 



9"K) 



sinh 



For small he, this relation simplifies to 
he = 



Ah 



(13) 



(14) 



and, for convenience, we shall use this expression for the 
proper distance to the brick wall. 



III. EXTENSION OF THE BRICK WALL 
MODEL TO HIGHER ORDERS IN THE WKB 
APPROXIMATION 

In this section, after a rapid sketch of the assumptions 
and approximations that are involved in evaluating the 
black hole entropy using the brick wall model, we go on 
to outline the procedure for computing the brick wall 
entropy at the higher orders in the WKB approximation. 



A. Basic assumptions 

There are two crucial assumptions in the brick wall ap- 
proach to black hole entropy. The first assumption con- 
cerns the modeling of the microscopic origin of the black 
hole entropy, and the second is regarding the handling of 
the divergences that arise close to the event horizon. 

As we have mentioned before, the brick wall model is a 
semi-classical approach wherein the black hole is assumed 
to be described by a fixed classical geometry. It is fur- 
ther assumed that the black hole is in equilibrium with 
a thermal bath of quantum matter fields at the Hawk- 
ing temperature. Moreover, it is the canonical entropy 
(actually, a specific component) of the quantum matter 
fields that are propagating outside the black hole horizon 
that is identified to be the entropy of the black hole. 

In the process of calculating the canonical entropy of 
a matter field outside the black hole horizon, we need to 
evaluate the density of states of the field. However, one 
finds that, due to the infinite blue shifting of the modes in 
the vicinity of the event horizon, the density of states ac- 
tually diverges. This divergence is regulated in the model 
by introducing a cut-off by hand above the horizon. The 
cut-ofi' — popularly referred to as the brick-wall — is basi- 
cally a static, spherical mirror at which the matter fields 
are assumed to satisfy, say, the Dirichlet boundary condi- 
tions. One finds that the leading component of the brick 
wall entropy diverges as h~^, where he is the proper dis- 
tance to the brick wall defined in Eq. (12). (The other 
component is essentially a volume dependent term that 
arises even in fiat space.) It is this contribution that 



is identified to be the entropy of the black hole. More- 
over, a specific choice for the cut-off he has to be made 
(this depends on the number of fields, the dimension of 
the spacetime, etc., but is generally of the order of the 
Planck length (.^^), in order to reproduce the Bekenstein- 
Hawking area law (2). As we mentioned, the area law (2) 
has been recovered in this approach for a variety of black 
hole spacetimes and matter fields [27 43]. 



B. Essential approximations 

Two approximations turn out to be essential to make 
the computation of the brick wall entropy tractable. The 
first approximation is required in evaluating the density 
of states of matter fields around black holes, and the sec- 
ond involves expanding the metric near the event horizon. 

As we pointed out above, in order to evaluate the brick 
wall entropy, one needs to evaluate the density of states 
of matter fields around black holes. However, apart from 
some lower dimensional cases, the density of states can- 
not be evaluated exactly. As a result, in the brick wall 
model, the density of states is usually evaluated at the 
leading order in h in the WKB approximation. 

Moreover, barring a few special cases, one finds that, 
even after the WKB approximation, the brick wall en- 
tropy cannot be evaluated exactly. Recall that the domi- 
nant contribution to the entropy arises due to the modes 
close to the horizon. Motivated by this feature, one Tay- 
lor expands the metric functions /(r) and g{r) near the 
horizon in order to obtain a closed form expression for 
the brick wall entropy. 



C. The methodology 

Having discussed the assumptions and approximations 
involved in the brick wall approach, in the remainder 
of this sub-section, we shall outline the procedure for 
evaluating the brick wall entropy at the higher orders in 
the WKB approximation. 

The key assumption of the brick wall model, as we have 
pointed out above, is that the black hole is in equilibrium 
with a bath of thermal radiation at the Hawking temper- 
ature of the hole. The free energy _F of a scalar field at 
the inverse temperature (3 is given by (see, for example, 
Ref. [12]) 



F = 







00 





dE (^^) ln[l-exp-(/Ji;)], 

T{E) 



(15) 



where T{E) denotes the total number of modes of the 

field with energy less than E. We have integrated the 
first of the above equation by parts to arrive at the second 
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and have assumed that the boundary term vanishes. The 
canonical entropy associated with the free energy F is 
given by 



and, it is this entropy, evaluated at the Hawking tem- 
perature, that will be identified to be the entropy of the 
black hole. 



(16) 



Consider a massive and minimally coupled scalar field $ that is propagating in the the line-element (4). Such a 
field satisfies the differential equation 



$ = 0, 



(17) 



where m denotes the mass of the field. The rotational symmetry of the line-element (4) allows us to decompose the 
normal modes Uj^^^, of the field $ as follows (see, for instance, Ref. [47]): 



R{r) 



j.D/2 Q\l1(j 



(18) 



where E, I and mj (with i £ [1, (D — 1)]) are the energy, angular momentum and the azimuthal angular momenta 
associated with the modes, respectively, the quantity G{r) is given by 



G{t) = V/W g{r) , 



(19) 



and Yimiid, 4>i) denote the hyper- spherical harmonics. On substituting the mode (18) in the equation of motion (17) 
and using the properties of the hyper-spherical harmonics, we find that the function R{r) satisfies the differential 
equation 



R"ir) + 



where the quantities V^{r) and A(r) are given by 



V\r) 
A(r) 



1 



G2(r) 
2G{r) 



- fir) 
AG^ir) 



-A(r) 



-I- 



D 



R{r) = 0, 



'{e + D-l)!f 



G'jr) 
G{r) 



+ 



D{D-2) 
4^2 



(20) 

(21a) 
(21b) 



The total number of modes T{E) of the field $ with energy less than E can be evaluated exactly if the solution to 

the differential equation (20) can be written down explicitly. However, apart from some simple (1 + l)-dimensional 
example[35], it proves to be difficult to obtain an exact analytical solution for the function R(r). As a result, the 
WKB approximation is almost always resorted to in the literature [27-42] , and it is the leading order WKB solution 
for R{r) that is utilized to evaluate the number of states T{E), and the resulting free energy F and the entropy of 
of the quantum field. Our goal here is to extend the analysis to the higher orders in the WKB approximation. 
Let us begin by expressing the function R{r) in the following WKB form: 



R{r) 



Co 



cxp 



r 

I J dfP{f) 



(22) 



where cq is a constant. On substituting this expression in Eq. (20), we find that the function P{r) satisfies the 
differential equation 



l\ (P"{r) 



A(r). 



^AJ \ P{r) J \2 J \ P{r) 
Let us now expand the function P{r) in a power series in h"^ as follows (see, for instance, Ref. [48]): 



P(r) = g?i2"P2„(r). 



(23) 



(24) 



n=0 
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On substituting this series in the differential equation (22) and collecting the terms of a given order in 7i^, we obtain 
following expressions for P2„(r) upto n = 3: 



Po(r-) = ±V{r) = ± 
P2(r) ^ ^ 

Mr) 
Pe{r) 



G{r) 

8Po{r)J \Po{r) 

2 



- fir) m2 + 



21 M 1/2 



4P^'(r)\ / A(r) 



5^*2'^^^ /^4P2(r)A(r) +P^'(r)>^ ^ /^SPa'W F'(r) - P2(r) ^"(r) 



2V{r) 



4F3(r) 



5P2(r)P4(r)\ /8P2^(r)+4P4(r)A(r)+Pi'(r)\ M(r)P2"( 



V{r) 



4F2(^) 



2F3(r 



2P^^(r)P2(r) +2P4(r) ^^(r) - 3 Pf (r) - 6 Pi(r) ^(r) 
8F3(r) 



(25a) 
(25b) 

(25c) 

(25d) 



Note that the function Pair) is related algebraically 
to the quantities Vir) and A(r). It is evident that the 
higher order functions P2nir) (with n > 0) can be ex- 
pressed in terms of the functions at the lower orders and 
their derivatives and, eventually, in terms of the func- 
tion Poir). 

On using the series expansion (24) in the standard 
semiclassical quantization procedure [12] , we can express 
the total number of states TiE) of the field with energy 
less than E as follows: 



TiE) = J2 ^2niE), 



(26) 



n=0 



where we have defined T2niE) as 



' fc2n-l - 

r2niE) = ( ) / dr lMi2i+D-l) 



xW(f)P2„(r), (27) 



with the quantity W(^) being given by 



Wil) = 



il + D-2)\ 
(£>-!)! l\ 



(28) 



It should be mentioned that, in the above expression for 

'^2niE), wc have approximated the sum over the angu- 
lar quantum numbers £ as an integral with a degeneracy 
factor W(^). Such an approximation is often made in 
the literature, and the approximation is considered to be 
valid since the separation between the states are expected 
to be small [33]. Moreover, the upper limit i^^^ on the 
integral over £ is a function of energy E of the mode 
and the radial coordinate r, and it has to be chosen such 
that Po(r) is real^. Furthermore, the lower limit on the 



integral over radial coordinate, viz. he, is the invariant 
thickness of the 'brick-wall' defined in (12), and the up- 
per limit L is the infra-red cutoff which we shall assume 
to be much larger than the horizon radius. 

A few clarifying remarks are in order at this stage of 
our discussion. In the semi-classical quantization of, say, 
a one-dimensional non-relativistic quantum particle, the 
integral over the coordinate will be carried out over the 
range wherein Pq is real [48]). In the case of bounded 
systems, these limits will prove to be the turning points 
of the potential, whereas in the case of potential barri- 
ers the limits will be between one of the turning points 
and infinity. In the context of black holes, the effective 
potential turns out to be a barrier and the integral over 
the radial coordinate is to carried out between the event 
horizon of the black hole (which is an infinity in terms of 
the tortoise coordinates) and the first turning point that 
is located on the barrier. But, one finds that, most of the 
contribution to the density of states of the quantum field 
arises due to the modes close to the event horizon of the 
black hole, while the upper limit located on the barrier 
leads to a volume dependent contribution to the entropy. 
As a result, the contribution to the number of states and 
the free energy and the entropy of the quantum field due 
to the upper (infra-red) limit is usually ignored in the 
literature. 

We should emphasize the point that, apart from re- 
placing the sum over £ by an integral, we have not made 

any approximations until now. Hereafter, wc shall make 
two approximations that we had discussed is some detail 
in the last subsection. Firstly, we shall approximate the 
line-elcmcnt (4) near the event horizon of a spherically 
symmetric black hole to be that of Rindler spacetime, 



^ Actually, the limits have to be chosen such that P2n{r) are real 



for all n. However, since, for n > 0, the functions P2n('') can 
be expressed in terms of Po(r) and the real functions V{r) and 
A(r), when Pa{r) is real, P2n(r) arc real as well. Therefore, the 
limits on t proves to be the same for all n. 
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viz. Eq. (10). It should be pointed out that such an ap- 
proximation is always made in the literature to arrive at 
closed form expressions for the free energy and the en- 
tropy of the quantum field. Secondly, we shall truncate 
the series (24) at a particular order (we shall work until 
the sixth order in h) , and evaluate the density of states 
and the associated free energy and the entropy of the 
quantum field around the black hole. It is important to 
note that, in the literature, it is only the leading term in 
the series (26) that has always been taken into account 
ignoring the higher orders when evaluating the brick- wall 
entropy**. 

In the following two sections, we shall evaluate the con- 
tributions to the brick wall entropy at the higher order for 
four and six dimensional black holes, respectively. As we 
shall see, the contributions to the entropy from the higher 
orders turn out to be of the same order as the leading or- 
der in the WKB approximation. In other words, it may 
be incomplete to calculate the contribution only from the 
leading term in the WKB expansion. Moreover, we show 
that the brick wall entropy leads to generic corrections 
to Bekenstein-Hawking entropy (2). For instance, in the 
case of four-dimensions, we find that the brick-wall en- 
tropy has the form: 

5r^ = ^BH+^^^°H^H)log(:^), (29) 

where J^^'^^^Ah) depends on the surface gravity and the 

second derivative of the metric at the horizon. 

Before we proceed with the calculations, there is yet 
another point concerning the WKB approximation at the 
sub- leading orders that we need to discuss. As we men- 
tioned above, the limits on the integral over £ has been 
chosen such that Pq (r) is real. This condition essentially 
identifies the turning points of the potential. Notice that, 
in Eq. (25), all the higher order WKB terms — i.e. P2n(?') 
for n > contain Po{r) in the denominator. Obvi- 
ously, these functions will diverge at the turning points, 
or equivalently, at the upper limits i. Such a divergence 
is a well-known feature of the WKB approximation at 
the higher orders [48], and we shall devise a systematic 
procedure to isolate these divergences. We shall outline 
this procedure in the next section and relegate some of 
the details to Appendix (B). 



IV. HIGHER ORDER CONTRIBUTIONS IN 
FOUR DIMENSIONS 

In this section, we shall evaluate the brick wall entropy 
for spherically symmetric, four dimensional black holes 



* The higher order WKB procedure we use is different compared 

to the approach used in the quasi-normal modes [49-51]. In 
Ref. [49], the Regge- Wheeler potential is expanded around the 
maxima, and the modes close to the maxima are matched to the 
one close to the horizon. 



by considering the contributions up to the n = 3 term in 
the series expansion (26) for the number of states of the 
quantum field. For simplicity, we shall consider here the 
case of /(r) = g{r) in the line-element (4) and restrict 
ourselves to a massless scalar field (i.e. m = 0). In 
Appendix (C), we shall extend the second order results 
we obtain in this section for the general case wherein 
/(r) ^ g{r) and, in Appendix (D), we extend the analysis 
to a massive field, but restrict ourselves to the case /(r) = 
9{r)- 



A. Second order 

Let us now evaluate the contribution due to the n = 1 

term in the scries (26). For /(r) = g{r), we find that the 
expression (25b) for second order 'momentum' P2 (r) can 
be written as 



+ A^(.) ( ^'^W 



(30) 



where the functions ('')) -^2 ('') P2 (r) are given 
by 



P^'\r) 



g'\r)\ /3g'(r)\ , /g"(r) 



g^r)"^ I \Arg{r)J \8g{r)J yAr"^/ ' 



+ 



y32j\g^{r)J \8rg{r)J \8r'^ , 
and, for convenience, we have defined 

g{£,r) = [£-X{r)f' (32) 
with £ = and A(r) being given by 

'airy 



(33) 



We now need to substitute the above expression for 
P2(r) in Eq. (27) and evaluate the number of modes r2 
with the upper limit ^max on the integral over £ being de- 
termined by the condition that the term g{£, r) vanishes. 
Clearly, the integral over £ will diverge in such a case. In 
order to isolate the finite contribution due to these higher 
order WKB modes, it is necessary that we follow a sys- 
tematic procedure. The procedure we shall adopt is as 
follows. We shall first rewrite all the terms containing 
inverse powers of Q{£,r) in terms of derivatives of £ as 
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follows: 



cussed in Appendix (B), we get 



1 



1 

1 



f dgi£,r) 



= -4 



d£^ 

3j \ d£3 



(34a) 
(34b) 

(34c) 



Then, before evaluating the i integral, we shall make use 
of the Leibnitz's rule, viz. 



b(x) 



a{x) 



= f[x,a{x) 



da{x) 
dx 



f[x,b{x)] 

b{x) 

dt 



db{x) 
dx 

df{x,t) 



dx 



(35) 



a{x) 



and interchange the order of differentiation and integra- 
tion over the energy E and £. When we do so, we find 
that the divergences occiir at the turning point. We have 
checked the procedure up to the 6'''^-order WKB modes 
and, indeed, systematically separates the non-divergent 
part from tlw^ divergent. For completeness, in Appendix 
(B), we give the details of the above procedure. (The pro- 
cedure involves calculating the contour integral around 
the branch cut that joins the turning points. For details, 
see Sec. (10.7) in Ref. [48].) 

Having obtained the non-divergent part of the mode- 
functions as a function of E, our next step is to evaluate 
the contribution of these modes to the density of states 
r2{E). Using the general expression (27), we have 



r2{E) = ^ J dr J d£{2i+l)P2{r). (36) 

r^+h 

Substituting for P2{r) from Eq. (30) and using the rela- 
tions (34), we get 



s 



1 f ,j'P^'\r) f^^dg{£,r) 



2 J " ■ d£ 

L e 



(37) 



^ J drr^P[^\r) J dX X 



d^g{£,r) 

d£^ 



1 3r2pf(r) 



+ - I dr - 'I J dX X^ 





d^g{£,i 

d£^ 



Using the Leibniz rule (35) and following the steps dis- 



r2{E) 



4rg'(r) 
Mr) 



+ r 



g'jrf 9"{r) 
3g{r)^ 2g{r) 



(38) 

Following points are worth noting regarding the above 
expression: (i) In the case of leading order WKB modes, 
the density of states goes as E^ [sec Eq. (A2)]. However, 
for the second-order WKB modes the density of states 
scales as E. (ii) As in the leading-order, most of the 
contributions to the entropy come close to the horizon, 
(iii) The expression for the density of state (C3) for the 
general spherically symmetric spacetime is same as for 
the special case discussed in this section. Hence, the 
dependence on the entropy with area is identical to the 
special case discussed in this section. 

Substituting the above expression in Eq. (15), and in- 
tegrating over E, the free-energy is 



dr 

+ h 



1 _ 4ryW 2 //Ml _ £>) 
3 3g{r) \3g{r)^ 2g{r) 



Using the relation (16), the entropy is given by 



(39) 



S2 = 



3hp 



1^ 

I 

„+h 



dr 



4r-ff'(r) 2 
3<7(r) 



g'jrf g"ir) 
3<7(r)2 2g{r) 



(40) 

As mentioned above, maximum contribution to the en- 
tropy is from the modes close to the horizon. Hence, 
using the expansion (11) close to the horizon and the 
definition of surface gravity (8), we get. 



IrJ, 



9"KK AC 

72 ^9" 



log ( # ) (41) 



where he is given by Eq. (14). This is the first result of 
this paper, regarding which we would like to stress the 
following points: 

1. The dependence of the entropy on area (from the 
second-order WKB modes) is similar to that from 
the zeroth order WKB modes (A6). Also the con- 
tribution to the entropy from the second order 
WKB modes contribute more as compared to the 
leading order WKB modes. This result has two 
immediate consequences: 

(a) To associate the brick-wall entropy to S^^ it 
is necessary to calculate all the higher order WKB 
mode contribution to the brick-wall entropy. 

(b) The sub-leading corrections (at the zeroth and 
second order WKB) depend only on the surface 
gravity and second derivative of the metric func- 
tions. They are of the form .F(^h) log(^H/^c)- 
To confirm the generic structure for higher-order, 
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in the next two subsections we evaluate fourth and 
sixth order contributions to the brick- waU entropy^. 

2. If the surface gravity is inversely proportional to 
horizon radius and g"{r^) is inversely proportional 
to the square of the horizon radius, then second 
term in the RHS of (41) is a constant. In this case, 
the corrections to S^^^ are purely logarithmic and 
docs not contain any power-law dependence. This 
uniquely corresponds to Schwarzschild spacetime. 

In the case of Schwarzschild, we have 



2M 

fir) = g{r) = 1 - — 



(42) 



where M is the mass of the black hole. The hori- 
zon is at = 2M, k = 1/(4M) and g"{r^) = 
— 1/{2M'^). Substituting the above expressions in 
Eq. (41), we get 



4M^ 



1 



(43) 



This result shows that, at least, in the zeroth and 
second order, there are no power-law corrections to 

S'gjj for the four-dimciisioiial Schwarzschild black 
hole, while, for all other black holes — since k and 
g"{r) has a more non- trivial structure - there are 
power-law corrections to the Bekcnstcin-Hawking 
entropy. This leads to the following conclusion: 
The power-law corrections to the entropy occur for 
any non-vacuum solutions. In Sec. (VI) we obtain 
the entropy for some known black hole solutions. 



B. Fourth order 

Using the expression (25c), we get, 
PM = (0 I Kr)Pi'\r) A^(r)pf (r) 
A3(r)Pf^(r) , \\r)Pt\r) 



+ 



G^{S,r) 



+ 



g^'{£,r) 



(44) 



where, the complete form of Pj*^ (r) [where i = • • • 4] 
are given in Appendix (E). 

Rewriting the above expressions in terms of the deriva- 
tives of energy and following the procedure discussed in 
Appendix (B), the contribution to the density of states 



^ It should be noted that, in the case of sixth order WKB modes, 

the integral over E is divergent near E ^ 0. However, the near- 
horizon contribution of the entropy in identical to the one ob- 
tained in this subsection. The fourth order WKB modes do not 
contribute to the brick-wall entropy. 



by the fourth-order WKB modes is given by: 



TiiE) = - I dr 

TT 





£ 

+ ^^i'^W^ / dXXgi£,r) 



£ 

- ^PPi^)-^ J d\\^g{£,r) (45) 







Integrating over A, we get 

„(4) 



r4E) = '-^ j drU^^Hr), 



(46) 



where Cq^'' is a constant and Y.'^^^r) is given in Eq. (E6). 
Using the expansion (11) close to the horizon, we get 



ME) = 



(4) 



323 r^K 5 rlg"{r)- 20 Kr^ 



2520 (r-rjj)2 



+ 



16{r-r^) 



(47) 

This is the second result of the paper, regarding which we 
would like to stress the following points: (i) The fourth order 
contributions to the density of states goes as 1/E. Using 
the expression (15), it is easy to see that the fourth order 
contribution to the free energy is independent of /3 and, hence, 
the contribution to the entropy vanishes®, (ii) The density of 
states contribution close to the horizon again depends only 
on the first and second order derivatives of the metric, (ii) 
Comparing the fourth order contribution to the density of 
states with the leading and second order, it is clear that the 
density of states scales as E^~^" where n is the order of the 
WKB modes. 



Sixth order 



Using the expression (25d), we get, 
P,^r) = I Hr)4'^(r) , A^(r)Pf (0 



e«(f,r) 



>?{r)Pr'(r) , A^(r)P6''^(r) A--'(r)P,^-^(r) 



G^^{e,r) 

>^\r)Pi''\r) 
g^'{S,r) 



Q^^{S,r) 



Q^^(S,r) 



(48) 



Note that, as mentioned earlier, the free-energy integral has an 
infra-red {E — » 0) divergence. 
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where, the complete form of Pg^^ (r) [where i = • • • 6] are 
given in Appendix (F). 

Rewriting the above expressions in terms of the derivatives 
of energy and following the procedure discussed in Appendix 
(B), the contribution to the density of states by the sixth- 
order WKB modes is given by: 



r,{E) = - 

TV 



dr 







10395 

256 
135135 

512 
^ 2027025 

Integrating over A, we get 










(49) 



6 l^Q^s J 



9 



(6) 7^ 

re{E)=^-^ J drE(«)(r) 

r„+h 



(50) 



where c^^ is a constant E'-®^(r) is given by Eq. (F8). Repeat- 
ing the steps i. e. using the relation (Eq. (15)) obtaining the 
free-energy, substituting the free-energy in (16) and expand- 
ing the metric close to horizon using Eq. (11), we get. 



Se _ 138 927r^ rl 
e ~ 45045 /i? 



977^,,. , 2 ^ 307r2 



(51) 



This is the third result of the paper, regarding which we would 
like to stress the following points: (i) The sub-leading correc- 
tions (like the zeroth and second order WKB) depend only 
on the surface gravity and second derivative of the metric 
functions. This indeed implies that brick-wall entropy does 
indeed provide generic corrections to the Bekcnstein-Hawking 
entropy at all orders. Wc have shown this to be the case upto 
sixth order. It is natural to expect this to be valid for all 
higher orders, (ii) As mentioned above, the density of states 
in each order is given by E^~". (iii) e in the above expression 
is due to the fact that the E divergences as £ — » 0. Thus, 
the above expression for the entropy depends on the infrar-red 
cutoff. 



V. HIGHER ORDER CONTRIBUTIONS IN SIX 
DIMENSIONS 

In this section, we obtain the zeroth and second-order 
WKB mode contributions to the brick-wall entropy in six- 
dimensional black hole spacetimc. The analysis can be ex- 
tended to any even dimensional spacetime, however, the anal- 
ysis in odd-dimensional spacetime is more involved'^. 

We show that the results of the brick-wall entropy in the 
zeroth and second order WKB modes have the same structure 
confirming the results of 4-dimensions and has the following 
generic form: 



4w' = s^^+g{Au)+:F^'''\Au) log 



(52) 



A. Zeroth order 

In the case of D = 4, the weight function (28) becomes 
{£+!){£ + 2) 



(53) 



Substituting the above expression in (27), the density of states 
for the zeroth-order WKB modes (for /(r) = g{r)) is 



■.(6D) 



L £ 



2 



where, g{£,r) is given by (32) and 



x = i{i + 3)re^. 



(55) 



Repeating the procedure discussed in the previous section, 
the zeroth order brick-wall entropy is given by 



c.(6D) _ 327r' 



L L 
5 /• 4 r „2 

dr——— + 



I 



Li 



945/95/15 ] g{r)^ 135/33R3 J g{' 

Expanding the metric near the horizon (11), upto third order, 
and using the relation (14), the zeroth-order entropy is given 
by 

= 37&+^-(-)+^-(-)l'^<|)'(^^) 

where 

^.2 

^<°) " 15120 fe2 ["^^"(''h) + le^r^ + 56] 
^(0) = +6^0 [(2«s"'(^h) - 3s"(r,)^) (58) 
-I- 24g"{r^)Krl - 224/t -|- {56g"{r^) - 48k^) r^] . 

As in the 4-dimensions, the leading order term in the above 
expression (57) is proportional to area (Bekenstein-Hawking 



^ This can be traced to the fact that the wave propagation in these 
spacetimes are non-local. For more discussion, see Refs. [52, 53]. 
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area relation). The sub-loading term has two parts: (a) One 
that contain a purely power-law corrections [5(o)('fH)] that is 
absent in the case of 4-(liuiensions. (b) The logarittnnic term 
contains a prefactor which, in general, is a function of area 
as like in 4-dimensions. J^^^j (^h), as in 4-dimensions, depend 
upto the second derivative of the metric close to the horizon 
while (rn) depend upto the third derivative of the metric 
close to the horizon. 



B. Second order 

Substituting Eq. (53) in Eq. (25b), the contribution to the 
density of states from the second-order WKB modes is given 
by 



r^+h L 


|pie^)(.)/dAW(A)A(.)^^ 



+ 4P«' 



(59) 



where, the degeneracy factor W(A) in terms of A is given by 

(60) 



mA) = ^^+2, 



and 



g{r) _ g'jr) 
3fl(r)2 3g{r)g'{r) , g'irf , g{r)g"ir) 



(61) 



+ ■ 



+ ■ 



4J.4 4^.3 8^.2 g^2 

5g{rf 5g(rfg'(r) , 5g(r)g'(rf 



8r6 



32'r 



Repeating the procedure discussed in Appendix (B), and sub- 
stituting the resultant in Eq. (15), we get, 



1^ 



+ 



— / 



dr 

TH+h 
L 

dr 



545(r)2 ^ 1085(r) 54 \ ' 



g"{ry + 6g'{ry r^g'{rf 2r^ 



1355(r)2 



90g(r)3 135t?(r) 



Substituting the above expression in Eq. (16) and expanding 
the metric using (11) the second order WKB mode contribu- 
tion to the brick-wall entropy is given by 



where 

6^(2) ('"h) 

•^(2) ('"h) 



'^^ + (rn) + :F,,, (rn) log ( ^ ) , (63) 



180M '^t^' ' 



^[80-r^5"(r,)] 



2160 

-i^^[{9"{r^f-2.g"'{r^))rl 



(64) 



-Sg"{r^)Krl - 800« - {112^ + 40g" {r^))r^] 



This is the fourth result of this paper, regarding which we 
would like to discuss the following: (i) As in the case of 
4-dimensions, the second order WKB modes contribute the 
same to the brick-wall entropy as the zeroth order modes. 
This again proves that in-order to associate the brick-wall en- 
tropy to the black hole entropy, it is necessary to calculate 
all the higher order WKB mode contribution, (ii) As in the 
case of 4-dimensions, the dependence of the entropy on the 
horizon area is the same in both the orders, (iii) Q^^^ (rn) (like 
^(0) (^h)) has a generic power-law corrections to S^^ and de- 
pend only upto the second derivative of the metric near the 
horizon. J^^^^{rn) (like ^(a^irn)) has a prefactor which is a 
function of the area, and — as in 4-dimensions - is a constant 
only for the Schwarzschild spacetime. 



VI. RESULTS FOR SPECIFIC BLACK HOLES 

In this section, we shall explicitly write down the brick 
wall entropy (evaluated upto the second order in the WKB 
approximation) for a few well-known black hole solutions in 
four and six spacetime dimensions. We shall restrict ourselves 
to the cases wherein f{r) = g{r). 



A. Four dimensional examples 

We find that, in four dimensions, on combining the zeroth 
order (A6) and the second order (41) terms, the total brick 
wall entropy can be expressed as 



(65) 



where, in order for the leading term to match the Bekenstein- 
Hawking entropy (2), we have set the brick wall invariant 
cutoff he to be 



QOtt 



(66) 



and the quantity T'^'^^^Ah) is given by 

^(-)(^h) = - (^) g"{T^) rl - (1) . (67) 



1. Schwarzschild black hole 

For the Schwarzschild black hole, the metric coefficients are 
given by Eq. (42) and the event-horizon of the black hole is 

located at t-j^ = (2Af). The surface gravity k and the second 
derivative of the metric at the horizon are given by 



1 



2A#2 



(68) 



On substituting these expressions in Eq. (65) , we obtain that 

(69) 
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2. Schwarzschild (anti-)de sitter spacetime 

For the Schwarzschild (anti-)de sitter spacetime, the metric 
function g{r) is given by 



air) = ( 1 - — ± 



r y 



(70) 



where y = {l/M)"^, r = (f/M), M is the mass of the 
black hole, I is related to the positive (negative) eosmolog- 
ical constant and — (+) corresponds to (anti-)de Sitter space- 
time. Note that the coordinates y and r are dimension- 
less. While the Schwarzschild anti-de Sitter spacetime has 
only one horizon associated with the singularity at the ori- 
gin, the Schwarzschild de Sitter has two — one event and one 
cosmological — horizons. Here, we shall focus on the entropy 
associated with the event-horizon. 

Recall that the event horizon is identified by the condition 
g{r) = 0. On substituting the resulting corresponding to 
the above g{r) in Eq. (65), we find that the brick wall entropy 
upto the second order can be expressed as 



,(4D) 

'Sch-(a)dS 



.1/2 



15 aT 



A 
ny 



log 



M'A„ 



(71) 



where A^ defined in-terrns of the coordinate r is also dinien- 
sionless. In contrast to the purely Schwarzschild case wherein 
the prefactor to the logarithmic correction was a constant, 
here the factor is a function of the horizon area. 



Reissner-Nordstrom black hole 



For the Reissner-Nordstrom black hole, we have 



g{r) = 1 



2M g; 



(r — r_)(r — r+) 



(72) 



where M and Q denotes the mass and the electric charge of 
the black hole. Also, r = r/M and r± is the outer/inner 
horizon given by 



1±A 1 



01 
M2 



(73) 



where, again, r is a dimensionless variable. It is the outer 
horizon r+ that is the event horizon of the black hole. 

On substituting the above relations in Eq. (65), we obtain 
the brick wall entropy upto the second order to be 



e(4D) _ „ 



.1/2 



15 Al/' 



log 



(74) 



where, again, A^^ defined in-tcrms of r is dimensionless. As 
in the previous example, the prefactor again turns out to be 
a function of the horizon area A^^ . 



where, as in the four dimensional case, we have chosen the 
invariant cutoff he to be such that the leading term matches 
the Bekenstein-Hawking entropy. The quantities he, G{Ah), 
J^'-'^°\Ah) are given by 



h„ 



:f^'°\Au) 



net 



(l8607r2 

(aSio) [(8'^^"'e^H)-5.-(..)).^ 

+ 40 g"{r^) KT'l + 2688 k 
+ 8 {46k^ +21g"{r^)) 



(76) 



(77) 



GiAu) = 



31 

165^4^ 



1/2 



V252/ 

X {-5 g"{r^)rl+ 8 Kr^+ 308) . (78) 



1. Schwarzschild black hole 

In six dimensions, the function g{r) for Schwarzschild black 
holes is given by 



9{r) 



(79) 



where r^ is related to the black hole mass (M) by the relation 



M ■■ 



(2jKrl 



(80) 



with Ge being the six-dimensional Newton's constant (which 
we shall hereafter set to unity). On using the definition (8) 
of the surface gravity n we find that 



3 

2r„ 



(81) 



where fjj = Vg . Substituting the derivatives of the above 
metric function g{r) in the expression Eq. (75), we obtain the 
brick wall entropy to be 



;,(6D) 



63 V 33 ei. 840 °^ I, £4, 




(82) 



where in deriving the above expression we have used the ex- 
pression for the area of the 4-sphere i.e. = {87i^)/3r^. 

Unlike four-dimension, there is a pure power-law correction 
term to the Bekenstein-Hawking entropy. 



2. Schwarzschild (Anti-)De sitter black hole 



B. Six dimensional examples 

On combining the zeroth order (57) and the second or- 
der (63) terms, we find that the brick wall entropy for six 
dimensional black holes can be expressed as 

4w ^ = 5b„ + 0{Au) + (.Ah) log (^^) , (75) 



The line-element for Schwarzschild (Anti-)De sitter space- 
time is given by (4) with 

/(r)=ff(r) = l-(^)'± J = l-1±- (83) 

where y,, = (l/vg)^, r {f/rfj, I is related to the posi- 
tive (negative) cosmological constant and — (-I-) corresponds 
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to asymptotic (Anti-)dc Sitter. Here again, for Schwarzschild- 
de Sitter, we consider only the event-horizon. 

The event horizon is given by the condition g{r = r^) =0. 
Substituting these in Eq. (75), we get. 



^(6D) 
'(A)dS 




1512 2/2 Wy, 252j/,rH 15rl 



Note that as in four-dimensions, is dimensionless. The 
above expression gives a series of power-law corrections to 



Reissner Nordstrom black hole 



For 6-dimensional Reissner Nordstrom black hole, we have 



where the charge of the black hole is given by 



Q 



3_e^ 

27r G„ 



r 



(85) 



(86) 



As in 4-dimensions, this has two horizons event and Cauchy 
horizon. The event horizon (r^) is the outer horizon while 
the inner horizon is the Cauchy horizon. Note that x and r 
are dimensionless. 

Substituting these in Eq. (75), we get, 



i(6D) 



•'^ 252 V 165 £1 



2 72 
308r„ + — 



234x 



(87) 



23x' X 
280ri2 I40r9 



3 - 28x 1 



840r6 



15r3 



log 



H 



where, again, A^ is dimensionless. As cam be seen, this also 
generates a scries of power-law corrections to the Bekenstein- 
Hawking entropy. 



where A = 6q^j^, r = ru; \ = \ lo and u> is related to 



the ADM mass (M 



ADM y by the relation 



— A 

4tt 



H 



(90) 



Note that the rescaled variables r, A are dimensionless. The 
horizon is given by the condition /(r) = and occurs at 



r = such that 



rl + 3ArH -1 = 0. 



(91) 



The existence of the horizon requires A > and which then 
gives rH^ < 1. The surface gravity of the event horizon is 
given by 



K = 



1/3 ri + 2rl 



2r„ 



2-r3 



(92) 



Substituting these in Eq. (75), we get, 

,,2/3, 



q(6D) 



125 



148955 w^'-V 

523908 ~ 
PI 

1 1175 



2hr± 



+ ■ 



31(2 ■ 



i) 



725 



31(2 
1 



216 (2 - 

355 1 
756 (2 - r 



r3 )6 



756 (2 - 

17 1 
189 2 -r 



,3)5 



504 (2 - 

655 1 



756 (2 - r3 )3 





649 


.3 14 
H' 


~ 7560 




\^—A 1 


log 









(93) 



It is instructive to stress the implications of this result: (a) 
In the limit of A ^ limit S^^^ reduces to S'^^^ . (ii) The 
sublcading corrections to the Bekenstcin-Hawking entropy for 
the Boulware-Deser black-hole have been obtained earlier by 
other authors (see, for instance, Ref. [58]). Using the Noether 
charge, it was shown that 



■^GB — 'S'bjj -I- 



87r2a;-^/3 



Ar„ 



(94) 



Comparing the two results we see that the subleading correc- 
tions in the brick-wall approach, unlike the Noether charge 
approach, can be completely specified by the horizon proper- 
ties. 



VII. DISCUSSION 



4- Boulware-Deser black hole 

Boulwaro-Dcser black-hole [56] is an exact spherically 
symmetric solution of the Einstein action modified by the 
quadratic Gauss-Bonnet combination i. e. 

1 = j d''x^ [r + (r' - mabR'''' + RabcdR"""'")] 

(88) 

where a^j^ is the Gauss-Bonnet coupling. The line-element is 
given by (4) where 



/(r) = g{r) = 




A. Summary 

As we have pointed out repeatedly, the brick wall model 
has been a very popular approach that has been utilized to 
recover the Bekenstein-Hawking entropy S-^^ in a multitude 
of situations [27-42]. In all these efforts, it is only the leading 
term in the WKB expansion (26) that has been taken into 
account in evaluating the density of states and the associated 
free energy and entropy of quantum fields around black holes. 
Also, the metric has almost always been assumed to be of the 
Rindler form near the event horizon. 

In this work, wc have extended the brick wall approach to 
the higher orders in the WKB approximation. Moreover, by 
expanding the metric functions /(r) and g{r) beyond the lead- 
ing order near the event horizon, we have been able to evaluate 
the corrections to the Bekenstein-Hawking entropy for spher- 
ically symmetric black holes in four and six dimensions. To 
begin with, wo have illustrated that, oven the often consid- 
ered zeroth order term in the WKB approximation leads to 
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corrections to the Bekenstein-Hawking entropy, provided the 
metric functions are expanded beyond the linear order near 
the horizon. Secondly, wc have shown that all the higher or- 
der terms in the WKB approximation have the same form as 
the zeroth order term. Lastly, we find that, the higher order 
WKB terms actually contribute more to the entropy than the 
lower order terms. 

Specifically, wc have shown that, upto the second order 
in the WKB approximation, the brick wall entropy of four 
dimensional black holes can be expressed as 

=^BH+^'^°'(-4H)log(^) , 

where T'^^'°\Aw) oc A"^ with n < 1. Whereas, in six dimen- 
sions, we find that the brick wall entropy up to the second 
order has the form 

SV = 5b„ +e(^H) +^'«°^(^h) log (^^) , 

where g{An) oc A"^ and J^('*°'(^h) oc A";^ with (n,m) < 1. 
Note that, while the brick wall entropy in four dimensions 
depends only on the first and the second derivatives of the 
metric at the horizon, in six dimensions, it depends on the 

third derivative as well. It is tempting to propose that, at 
least in even dimensions, the brick wall entropy will depend 
on as many as derivatives of the metric as half the number 
of spacetime dimensions! However, the black hole entropy 
is a coordinate invariant concept. If the brick wall entropy 
depends on arbitrary derivatives of the metric functions at 
the horizon, then it is not apriori evident that the resulting 
entropy will be coordinate invariant. Wc believe that this is 
an issue that needs to be addressed satisfactorily. 



B. Comparison with results from other approaches 

Power law and logarithmic corrections to the Bekenstein- 
Hawking entropy Sgn that we have obtained in the brick 
wall approach has been encountered earlier in a few other 
approaches to black hole entropy. For instance, the Noether 
charge approach predicts a generic power law correction to 
the Bekenstein-Hawking entropy [15]. However, unlike our 
approach wherein the brick wall entropy can be completely 
expressed in terms of the metric and its first few derivatives 
at the event horizon, the Noether charge entropy can not be 
mapped to the horizon properties. It is also interesting to note 
that, in the case of the four dimensional Reissner-Nordstrom 
black hole, for large horizon area, i. e. when M ^ £pj, the 

brick wall entropy S'^jJ^ [cf. Eq. (74)] reduces to 

5rn -Sbh- (^^J (^-^-^j^J • (95) 

Similar power law corrections arise on evaluating the entangle- 
ment entropy of such black holes [24] . This behavior seem to 
suggest a possible relationship between the brick wall model 
and the approach due to entanglement entropy. Another in- 
teresting feature is the absence of power law corrections in 
case of four dimensional Schwarzschild black hole. It seems 
to indicate that power law corrections to the Bekenstein- 
Hawking entropy arc related with the presence of matter. The 
logarithmic corrections that we have obtained as in Eq. (69) 



for the case of the four dimensional Schwarzschild black hole 
has also been arrived at in other methods such as the approach 
through conforrrial field theory [19], statistical fluctuations 
around thermal equilibrium [23] and spin foam models [18]. 
However, it should be pointed out that the prefactor to the 
logarithmic term that we obtain turns out to be different from 
the one that arises in the other approaches. 



C. A few words on the divergences 

The divergence that results in the need of a brick wall cut- 
off arises even at the leading order in the WKB approxima- 
tion, and is, obviously, well-known. So, it is not all supris- 
ing that such a divergence occurs at the higher orders terms 
as well. However, in addition to the brick wall divergence, 
we seem to encounter three more types of divergences at the 
higher orders. The first is the divergence that occurs at the 
upper limit when integrating over t in the higher orders and 
the second is an infra-red divergence that arises at a suffi- 
ciently high order when integrating over E (as in the case of 
the sixth order term in four dimensions). As we mentioned 
above, the higher order terms contribute more to the brick 
wall entropy than the lower orders. Therefore, the third and 
last is a divergence that can arise if we get around to summing 
over all the terms n. 

The first of these additional divergences is associated with 
the turning points. Such divergences arc known to occur at 
the higher orders in the WKB approximation even in non- 
rclativistic quantum mechanics. Evidently, these divergences 
are not a feature of the field theory, but a feature of the 
approximation. The procedure we have adopted to isolate 
and discard these divergences effectively deals with them. 

In contrast, the remaining two divergences that occur at 
the higher orders are field theoretic divergences. The infra-red 
divergence is clearly one as it arises when integrating over all 
the modes. In the context of the leading order results, it has 
been argued that the brick wall divergence can be absorbed 
into the renormalization of the Newton's constant [30, 41]. 
Clearly, if the higher order WKB terms continue to contribute 
more to the brick wall entropy than the lower order ones, then 
a divergence will arise when the contributions from all the 
orders are summed over. We believe that, when working at 
the higher orders in the WKB approximation, these additional 
divergences need to be accommodated in a renormalization 
procedure, along with the brick wall divergence itself. 



D. Outlook 

Since the odd dimensional cases are analytically more in- 
volved, after first working in four dimensions, we had jumped 
to consider six-dimensional black holes. Needless to add, it 
will be interesting to extend the current analysis to black holes 
in odd dimensional spacctimcs. The Banados-Teitclboim- 
Zanclli black hole in three dimensions [54, 55] and the five- 
dimensional Boulwarc-Deser black hole [56] are interesting 
cases that are to be studied. The canonical entropy has been 
calculated exactly around the the Banados-Teitelboim-Zanelli 
black hole (see, for example, Ref. [57]), and the entropy of the 
Boulware-Deser black hole is expected to contain a power-law 
correction (cx A\(^) to the Bekenstein-Hawking entropy (see, 
for instance, Ref. [58]). It will be worthwhile to investigate 
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as to how the brick wall entropy compares with these results. 
We hope to consider these cases in a future publication. 



However, if we expand the metric to higher orders (11), we 
get 
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APPENDIX A: 



RECOVERING S^j, FROM 



BRICK-WALL 

In the first part of the appendix, we provide the steps lead- 
ing to the Sgjj in the Schwaxzschild-like coordinate system i. 
e. /(r) — g{r) in the line-element (4). In the second part, we 
obtain the same in the tortoise coordinate (5). As mentioned 
in Sec. (Ill), the brick- wall model assumes that the the WKB 
mode functions are a good approximation for the radial modes 
near the horizon. In the case of Schwarzschild-like coordinate 
system, it is not apparent whether such an approximation is 
valid. 



1. Schwarzschild-like coordinate 

In the case of massless scalar field, the leading order WKB 
modes are given by 



9{r)- 



1/2 



Substituting the above expression in Eq. (27), we get 

L 



r dr 



(Al) 



(A2) 



Substituting the above expression in (15) and integrating over 
E, the free energy F now reads 



Fo 



and the entropy is 



2n' 



1^ 



r„+h 



dr , 



Stt-' 1 

45^^ j g'^ir) 

r„+h 



dr . 



(A3) 



(A4) 



On expanding the metric near the horizon up to the first- 
order, we recover the following standard result [12] 



So. 



90/i? 



+ 



90 



360 



log 



hi 



(A6) 



It is interesting to note that the form of the brick-wall entropy 
in the lowest-order WKB is same as the one obtained in the 
higher-order WKB. See Sec. (IV) for details. 



2. Tortoise coordinates 

For the case /(r) = g{r) in (4), the tortoise coordinate 
simplifies to 



-I 



dr 
9{r) 



(A7) 



As mentioned in Sec. (II), in the Rindler approxima- 
tion (10) the new radial coordinate x takes the form x = 
\og(r — r^)/g'{r^). Thus, the horizon r = corresponds 
to a; — > — cx). Hence the tortoise coordinate cover the range 
(—00,00). The massless scalar field propagating in the back- 
ground (5) is nothing but the Regge- Wheeler equation given 
by 



d^R 



-^+[oj'-VnwHx)\]R = 



(A8) 



where 



VR.w[r(x-)] 



The Regge- Wheeler potential decays exponentially near the 
event-horizon and as a power-law near spatial infinity (for 
asymptotically fiat spacetimes), i.e.. 



V[r{x)] ' 



exp{2Khx) ]V[r{x)] 



V — >-oo 1 



(AlO) 



Thus, the general solution to Eq. (A8) as a; — » =Foo can be 
written as a superposition of plane-waves: 



R[x\ &q){iu)x) -h C2 exp(— iwa;) , 



(All) 



where Cf^Ct are the constants determined by the choice of 
the boundary conditions. Thus, in other words, the WKB 
modes is a good approximation for the radial modes close to 
the horizon. 

Using the procedure discussed in Sec. (HI), the leading 
order density of states is given by 



r(->(£) 



37rR3 



1^ 



9{r)g'{r) 



-| 3/2 

, (A12) 



where H 'lsa, large positive number corresponding to the cutoff 
h. At the linear order in the near horizon approximation we 
finally obtain. 



/ 



dx 



g'K) 



^ exp {-2kx) (A13) 
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Following the procedure discussed above, the leading order 
free energy and the entropy are given by 



3 2 

exp {2kH) 



^«exp(2K/f). 



(A14) 



(A15) 



Now, we want to express this entropy in terms of invariant 
cutoff He defined as, 



Hc= J dr* y/g{r) = y ~ ^^P . 



(A16) 



Using this invariant cutoff, the final expression of entropy is 
given by, 



_ ' h 

° ~ 9om 



(A17) 



The entropy obtained in the tortoise coordinate is same as 
obtained in the Schwarzschild coordinate. 



APPENDIX B: ISOLATING THE FINITE 
CONTRIBUTION 

In this appendix, wc shall outline as to how we isolate the 
finite part of the integrals using the Leibniz rule (35). 

The first integral in the RHS of Eq. (38) does not lead to 
any divergent term. Using the relation (35) — with a{£) = 

0, b{£) = £ — we get 



I 



dg{£,r) 
d£ 



I d\P^°\r)G{£,r). (Bl) 



The second and third integral in the RHS of Eq. (38) lead 
to divergent terms. The origin of the divergent terms can be 



associated to the break-down of WKB approximation at the 
turning point {£). In order to see this, we evaluate the two 
integrals explicitly. 

Using the relation (35) in the second integral of Eq. (38) , 
we get, 



I 



d£^ 



(B2) 



d f\ dgi£,r) , dg{£,T) 



d£^J 



\{r)g(£,r)d\- 



£ 



(f - A)i/2 



where we have used Eq. (Bl) to obtain the second equation. 
Using the relation (35) in the third integral of Eq. (38), we 



get 



/ 



,2 d^g{£,r) 
d£^ 



d\ = 



d£3 

d_ 

d£ 



jd\\^g{£,r) (B3) 



£^ 



2e(f,r) 



£^ 



4gH£,r) 



From Eqs. (B2, B3), it is clear that both the integrals 
have a finite and divergent part. The divergence occurs at 
the turning point £ = X. This is not a physical divergences, 
this is occurring due to the fact that the WKB approximation 
is not valid close to the turning points. However, it can be 
shown that by introducing a cutoff close to the turning point 
the results arc independent of the cutoff. (For details, see Sec. 
(10.7) in Ref. [48].) 



APPENDIX C: SECOND ORDER CONTRIBUTION WHEN f (r) ^ g(r) 

Earlier, in Subsection IV A, we had computed the brick wall entropy at the second order in the WKB approximation for the 
specific case wherein /(r) — g{r) in the line-element (4). In this appendix, we shall obtain the corresponding result for the 

more general case of /(r) 7^ g{r)- For simplicity, wc shall again consider a masslcss field. 

When /(r) 7^ <?(»■), we find that, the 'momentum' at the second order P2{r) [cf. Eq. (25b)] can be written as 

P,»V,^(t^).M,.)(if^)..>,(|M), 

where g{£,r) is given by Eq. (32). We have defined the functions Pzoir), -Picl'') and P^'§{r) to be 



G(r)/'(r)\ fG(r)g'{r) 



^^(^''^^ 4r4 ) [ 8r-' ) ^ \ 16 r-\f{r) J \ 8r^g{r) J ' V 16r^g{r) J ' y 8r^ 



4r/(-r) / \ 4rg(r] 

:iG(r)/(r)\ /5G(r)/'(r)\ , fG{r)f'irf\ f G{r) f{r) f/ (r)\ , f G{r) g' ir)\ , /G(r)/"(r) 



8r« J \ 8r5 / V 32r4 

with the quantity G(r) given by Eq. (19). On we find that the number of states at the second order r2'^\E) can be expressed 
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as follows: 



2 3 \ 



r„+h 



(C3) 

Repeating the procedure discussed in the last appendix to identify and ignore the divergences in the above expression and 
substituting this expression in the integral (15) for the free energy, we obtain that 



-^2 — 



C 



-3rV(r)/'(r)5'(r)-6rVMfl(r)/"(r)] . (C4) 

Using the relation (16) and expanding the functions /(r) an g{r) about the event horizon as in Eq. (11), we obtain the the 
expression for the entropy to be 



9"{rn)\ , ff"{rJg'K)\ , f^^r^\ fl'(r„) 



144 



+ 



48 /'(r„) 



+ 



9pn J V /'(t-h) 



(C5) 



APPENDIX D: SECOND ORDER CONTRIBUTION FOR A MASSIVE FIELD 



In Subsection IV A, we had evaluated the brick wall entropy at the second order for a massless field in four spacetime 
dimensions. We had considered the specific case wherein /(r) = g{r) in the line-element (4). In this appendix, we shall discuss 
the corresponding result for a massive field. 

For a massive field, the quantity P^ir) as given by Eq. (25b) can be expressed as follows: 



Ptir)\ /^A(r)pW(r) + Pi^(r)\ ( A^W Pj^^r) + A(r) [P^ (r) + Pj^^ (r)] ' 



QmH£,r) 



(Dl) 



where we have defined the functions Pim(^)' P2m('') ^'^'^ ^imC'') to be 
^ - ^ ■ ''^m? g{r)g'{r) 



P^^r) = [g'\r)+g{r)g"{r)\ , P«(r) = 

with the quantity Qm{£,r) being given by 



16^3 j ,^g{r) + rg'{r)\ and Pi^(r) = (^^j [m' g{r) g'\r)\ , 

(D2) 



Gm{£, r)=[£ - A(r) - rn^gir)] 



1/2 



(D3) 



Carrying out the procedure we had discussed in Appendix B to identify and discard the divergences, we obtain the second 
order density of states rj™' {E) for the massive field to be 



+ 4Pi:2M ^ jd^ G(r.){£,r) + Pf (r) ^ J dX X' g^^)i£,r) 




I) [P^2 (r) + P^'2 (r)] ^ jdXX Q^^) (f, r) 



(D4) 



where £m =\f- ~ fl'(^)] • As in the massless case, on expanding the function g{f) near the event horizon as in Eq. (11), we 
find thatr^"'^(S) can be expressed as 



(D5) 



18 



Note that the last term containing the mass m of the field is inversely proportional to E. Recall that, in Subsection IV B, we 
had encountered such a behavior at the fourth order for the rnasslcss field [cf. Eq. (47)]. As we had then pointed out, such 
a dependence on the energy E in the number of states leads to a free energy that turns out to be independent of the inverse 
temperature (3 and, hence, the term does not contribute to the entropy. Therefore, the massless and the mgissive fields lead to 
the same entropy at the second order in the WKB approximation. 



APPENDIX E: EXPLICIT FORMS OF Pi'^(r) 
The functions Pj'^ (r) [where i = 0, • • • ,4] are given by, 

-^5(r)ff'(r)P^(°)(r) - \g{r)P^°\r)g"{r) - iffM^P^^W , (El) 

Pi'\r) = -5gir)P^°\r)P^'\r) + ^ [girfg'ir) P^'\r) + g{rfp;,^'\r)] 

[9{r)P^°\r)g'{rf+g{rfg'{r)P^^°\r)\-^g{r)g'{r)P^'\r) 
-\ [9'{rfP^'\r) + 3g{r)g'ir)Pi^'\r)+g{r)g"ir)P^'\r)+girfP^'^'\r)] , (E2) 



Pf )(r) = -5g{r)P^'>ir)P.^\r) - ^g{r)P^'\rf - ^g(r)g' (r)P.:^'\r) 



-\ [ff(r)s"(r)Pi'^(r)+s(r)V;(^\r)+s'(r)^Pi^^r)] - ^g{rfp^''\r) 
+ 4^ffWV(r)Pi°^(r) - ^ ^g{TfP^^\r) - \g{Tf g' {rf Pf\r) 

^3 [l5<?(r)V(r)Pi'^(r) +93(r)^^Pf )(r)] - \g(j-)g\T)Pf\r) 
^ [llff'(r)^5(r)Pi^'(r) + 9ff(r)^ff'(r)Pf ^r) + 2g{rf ^' {T)P^^\r)\ (E3) 



+ 4^.3 



Pf\r) = -5gir)P^^\r)PPir) - g,.(r) V^^) (.) + 23,(.)3p(^) (.),'(.) 



4 UgirfP^'Hr) - ^g{rfg'{rfP^'Hr) 



16^ 



+^ [25g{rfP^'\r)g'{r) + 13g{rfP^^'\r)] 

[l7g{r)PPir)g'{rf + 13g{rfg'{r)Pi^^\r) + 4ff(r)^Pf )(r)ff"(r)] , (E4) 



Pt\r) = -^girrP^\r) + ^,(r)V(r)Pf '(r) - lg{r)P^\rf - ^girfg'irfp^\r) , (E5) 



E(4)M = 323 r^g'jrf lOlr' g' {rf g" (r) - 631 r g' {rf 
^ ' 10080 3(r)2 1680 5(r) 

+S V^^^^ + V(r)ff'''(r) + 5sW(r)a(r)] + ^ [ihhg' {r)g" {r) + 2h2g^''\r)g{r)\ 

467g'(r)" + 150ff" (r)g(r) 17 g(r)^ 1223 ff'(r)ff(r) 

420 630 r2 2520 r ^ ' 
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APPENDIX F: EXPLICIT FORMS OF P^''(r) 
The functions Pg'' (r) [where i = 0, • • • ,6] are given by, 

P^M = -25(r)Vr(r)^-55(r)Pr(r)Pr(r)-^(!|>l [ffMPf 

+2g{r)g"{r)P^''\rf + 2g{r)P':^''\r) + 2g" {r)Pi''\r)\ - \pi''\r)g' {rf (Fl) 



^e^W = -^fW [^M^-rM' + ZPrw] [ff(r)fl'(r)Pf (r)Pi^)(r)+fl'(r)Pii)(r)] 

[25(r)Pr (r)^ff'(r) - ISp'^Pf (r) + ff(r)^Pr WP^^r) - 9ff(r)Pr (r)] 
+1^ [ff(r-)V(r-)J'r«J^2*°\0 - ll5'«'Pi°^(r) - 9g{r)g' {r)PT\r) - g{rf g" {r)P^°\rf] 
-^9{rf9"{r)PT\r) - Qg{rf P^°\rf Pi'^\r) - hg{r)Pi^\r)Pi°\r) - 5ff(r)Pf (r)Pf ^(r) 

+ '-^[9'(,rfPr\r)P^^\r)+Zg{rfPr(r)P:}''{r^ 

-2g{r)g"{r)P^'\r)P^'\r)\ - \ [g\rfPi'\r) + girfP^'-'^r) + g{r)Pi'\r)g'\r) + Mr)9\r)P',^'\r)\ (F2) 



Pf>(r) 



[3,MPr + 23Pr m] + 
[.MPrWPi^\.) + pf\r)] - 

[pf)w + 2s(r)Pr(r)Pf (r)+sWs''(r)P«(r)^] [SP^M 

\2 



[3ff(r)Pi">(r)2+46Pf'(r)] 

[35(r)Pi°\r)V46Pr(r)] 



+ 



4r 



^ [85(r)5'(r)Pi°^(r)Pi^>(r) + 25fl'(r)Pi^> (r) - 7g{rf P^'\r)P^'-''\r) 



+ [-175'(r)Pi^^(r) + 7g{rf P^'\r)Pr {r) - MrfP^°\r)P!,''\r) - 13s(r)Pl(^)(r)] 

_i(rmi [Pf WPi^\.) + P,«(.)] - [9{r)P^'\rf + 2gir)P^^\r)P^^\r) + Pfw] 

-6g{r)' [P^°\r)P^'\r)' + PfW^f '(r)] - 5g{r) [p^''\r)Pi°\r) + P^'\r)Pi'\r) + P^°\r)Pi^'\r)] 

9ir)9'(rr [p(i)(^)2 ^ p(o)(,)p(2)(,) _ -^pf (r) - -l^Pf >(, 



(F3) 



P^°\r)P^''-^>{r) + P^'>{r)P^'''>{r) + 



y/(2) 



(1)/- 



d"(2) 



M + P^'Hr)P^'^'\r) + lp;^'\rf + ZP'^ {r)P',^'\r) 
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P«(r)^+2Pr(r)Pf(r) + ^Pf)(r) 



ff(r)'ff'(r) 



13Pf \r)Pr'(r) - lPr{r)Pr\r)^Zg{rTP^'\r)P'^^'\r) - 



y(o) ^ 



j(0) 



'(2)/ 



4^3 

ffWV(r) 

4j,3 



^4p(l)^ 



16r4 



y(i)/ 



[95(r)Pf ^ (r)Pi^' (r) + 49Pi^' (r)] 
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y(2) 



(r) 



^Pi^'(r) + 6Pi^'(r)^ + 12Pr (r)Pf ^(r)' 



[3<;'(r)Pf)(r)+ff"(r)Pf (r)] 



Pi'-\rf + 2pW(r)Pf )(r) - ^^^f ^(0 



^(^[23,'WPf)W + 17,WPf^W] 

[i3Pf ) (.)Pf ) (r) + 3P,« (r )Pf ^) (r) - 7Pi"' (r)Pf ) (r)] 
[s(r)Pi^^(r)Pi'\r) + Pf ^(r)] - 2g(r)^Pi^^(r)^ - 125(r)2pf (r)Pi^'(r)Pi')(r) - is(r)^Pf ^'(r) 



8r2 
9(r)g'(r) 



-5g{r) [p^'\r)Pi'\r) + P^'\r)Pi''\r) + P^''\r)Pi'\r)] 



+ 



9ir)g\r) 



[g{r)Pi''\r)Pi^\r)-Pi^\r)\ 



+ ^ [3P^(i)(,)P^(2)(,) _ 2Pii)(r)P;(^'(r) - 2Pi^)(r)P;(i)(r) - P^'\r)P'^'-'\r)] 



(F4) 



Pi^^(r) = - 



8r6 



9{rfg"{r) 



+ 



g'(r)g(r)^ 

g(r)^g'(r)^ 
32r4 

4 J. 3 
g(r)V(r) 



46Pi°^(r)Pr (r) +3P2^^'(r)^ + ^Pr(r) 



3(2) 



(1)/ 



n2 , 166 ^(2) 



gW 



63(r) 



3P«(rf +46Pr(r)Pf(r) + ipP, 



(2) 



(r) 



9'{r)g{r) 
2r 



[9{r)P^'\r)P^^\r) + Pi'\r)] 
g{r)P^''\r)P^'\r) + Pi'\r)] 
[gir)P^^\rf + 2Pt\r)] 



[mg{r)g'{r)P^'\r)P^'\r) + 45g{r)P^'Hr) + 845(r)Pi(^^' (r)] 
[9Pi^'(r)P^(^'(r) +s(r)V,W(r)P^(^'(r)] - [g(r)Pf >(r) 



^+Pr 



9Pf )(r-)Pf - Pi^'(r)P^(^)(r) - J^Pf ^(r) - 29^Pf (r) 



8^2 L 2 V/'2 -2 V/-2 

-65(r)^ [pi^nr)'Pi''(r) + P^°\r)P^''Hrf] - 5gir) [p^''\r)PtHr) + P^'Hr)Pi'\r) + P^°\r)Pt\r)] 



+ ^g(r)g'(r) \g'{r)P^'\r) 



2^P 
9{r) 



(4) 



(r)-6P;<"'(r 



+ l9{rf 



P^^''\rf-2P^^>{r)P^'^\r 



3"(4) 



gW 



(F5) 
(r) 



Pf'(r) = 



[5(r)Pi^)(r)Pf)(r) + 25Pf^(r) 



_5g(ry 
4^6 

-^g(rf [gWPi'^WV2Pf)( 
5g(r)'3'(r) 



5g'(r)ff(r)^ 
4^.5 



[g(r)Pii'(r)Pf)(r) + 25Pf' (r)] 



16^- 



|lg(0'g'(0' [g(r)Pi''(r)Pi''(r) 



+ 25Pf' (r 



4r3 

35g(ryM!p(4) 

8r2 



2g{r)P^'\rf + llgirfPtHr)] - ^ [g{r)P^'\r)Pr\r) +25Pi^'\r)] 

[9{r)P^''\r)P^^^\r) - 5Pi^*\r)] - Qg{rf P^^\r)P^^\rf 



Pr(r) + 
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8r2 



-^^^^^f^ [9{r)P^'\rf - 2Pt\r)] - 5,(r) [P<^'(r)Pf '(r) + P^^'MPi^'w] 



(F6) 



Pi«^(r) 
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8r6 



g«'[g«Pf^M'-70P«M] - 



5g(r)-^ff'(r) 



[g{r)P^^\rf -7QPi^\r)] 
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32r4 



[5(r)Pf)(r 



70P 



(4) 



(r)] - 2girfP^^\rf - 5fl(r)p(^'(r)Pf '(r) 



(F7) 
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= 9341g(r)" A7Alg\r)g{rf 1308784 g"(r)g(r)« + 3926504 g'(r)^g(r)^ 
^ 180180 r4 16380 r3 5765760 

536120g<^) (r) g{rf + 3719032 g'(r) g"{r)g(rf + 2869040 g'jrfgjr) 

5765760 r 

, 213928g'^) (r)g(r)^ + 761280g" (r)^g(r)^ + 1261208g'(r)g<^>(Ofl(^)^ + 1508748g' (r)^g"(r)g(r) + 435674g'(r)'^ 

5765760 

' fl37280gf^H^)ff(^)^ + 947804g" (r)g(^'(r)g(r)^ - 

+ - 



r (l37280g''^^r)g{rf + 947804g" {r)g^^\r)g{rf + 903188g'(r)g(*)(r)g(r)^ + 462228g' (r)g"(r)2g(r)) 



+■ 



5765760 

r (971568g' {rfg'-^\r)g{rf - 4496g'(r) V'(Off(^) - 24496g'(r)^j 
5765760 g(r) 

(3473g'{rf - 27895g(r)g"(r)g'(r)^ + 84032g(r)^g'^^(r)g'(r)^ + 113100g(r)^g"(r)^g'(r)' 

5765760g(r)2 

(l83898g(r)3gW(r)g'(r)2 + 316316g(r)3g"(r)g(3)(r)g'(r) + 99528g(r)^g(5)(r)g'(r)) 

5765760g(r)2 

(40040g(r)3g"(r)3 + 61776g(r)*g(3)(r)2 + 123552g(r)*g"(r)gW(r) + 12012g(r)Sg(«)(r)) 

5765760g(r)2 
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